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A family of unitary Q-Ensembles of random matrices with governable confinement potential 
V (x) ~ I a; I" is studied employing exact results of the theory of non-classical orthogonal polynomials. 
' The density of levels, two-point kernel, locally rescaled two-level cluster function and smoothed 

, connected correlations between the density of eigenvalues are calculated for strong (a > 1) and 

border (a = 1) level confinement. It is shown that the density of states is a smooth function for 
Q > 1, and has a well-pronounced peak at the band center for a < 1. The case of border level 
O , confinement associated with transition point a = 1 is reduced to the exactly solvable PoUaczek 

' random-matrix ensemble. Unlike the density of states, all the two-point correlators remain (after 

proper rescaling) to be universal down to and including a = 1. 
PACS number(s): 05.45. -l-b, 05.40. 
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I. INTRODUCTION 

o . 

' Random matrix theory ||] pretends to description of a great variety of physical systems from complex nuclei and 
classically chaotic systems to electron transport in disordered conductors. The latters are known to exhibit crossover 
■ from metallic phase with Wigner-Dyson level statistics to insulating one with Poisson level statistics as disorder 
increases, passing through the Anderson transition where a third universal statistics should exist 
"Y*^ ■ There are two ways to describe such a crossover within the random matrix theory. The first one deals with random 
' matrix ensembles which do not hold their invariance under orthogonal, unitary or symplectic transformations from 
, the beginning due to the presence of a symmetry-breaking term in the joint probability density function P [H] of 
I ' the N X N elements of the random matrix H. Such random matrix ensembles with primordially broken symmetry 
' ^ ] demonstrate a deviation of the level statistics from the Wigner-Dyson type to the Poissonian one. 

The second approach starts with invariant ensembles of random matrices that implies the logarithmic repulsion 
between the levels confined by a parameter-dependent potential. In this case [|l| a rigorous treatment in the terms 
of orthogonal polynomials turns out to be successful if those are known for the chosen confinement potential. For 
instance, in the model based on the non-classical q-polynomials Q| the relevant parameter, entering the confinement 
potential, is associated with the strength of disorder. The strong confinement was found to be relevant to the weak 
$H ' disorder (metallic regime), while soft confinement potential was characteristic for strong disorder (insulating regime). 

Correspondence between these two approaches has been traced in where it was shown that in random-matrix 
ensemble with quadratic logarithmic confinement potential the spontaneous breaking of underlying symmetry of P [H] 
occurs. The symmetry breaking manifests itself in the loosing of the translational invariance of the two-level cluster 
function I2 (s, s') and leads to the Poisson-like level statistics. 

Another family of unitary ensembles of random matrices with governable confinement potential was proposed in , 
where so-called a-Ensemble was introduced. Symmetric a- Ensemble was considered in This ensemble provides an 
excellent basis for rigorous treatment Q and enables to explore how soft the confinement potential must be to cause 
deviations from the Wigner-Dyson statistics. 

Let us consider a physical system with broken time-reversal symmetry described by an A'^ x A'^ random matrix H 
with eigenvalues {x„} ,n — 1, A^. The joint probability density function P ({x}) can be written in the form |l| 

P {{x}) = Z-^ exp (^-(3 (^V {x,) -J^^\n\x, ~ x,\]i^ , /3 = 2, (1) 

which implies a pairwise logarithmic repulsion between the levels confined by the potential V (x) ; Z is a partition 
function. 

Symmetric a-Ensemble is characterized by parameter-dependent potential 
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Vc {x) = \ \x\ 
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supported on the whole real axis x G ]— oo, +cxd[. 

The Monte Carlo simulations 0, based on the mapping of initial random- matrix problem onto the problem of 
interacting particles in confined ID plasma, showed that in such a-Ensemble the deviations from the Wigner-Dyson 
statistics may also occur, and they take place near the center of the spectrum for < a < 1. The value a = 1 was 
associated with the sharp transition point corresponding to the crossover from strong to soft level confinement. Just 
at this point, a = 1, the mean- field approximation widely used in the random-matrix theory fails, giving rise to the 
singularity of the level density near the spectrum origin. 

In this paper we present the analytical treatment of a-Ensemble. In Section II a-Ensemble with strong confinement 
potential (a > 1) is studied employing exact results of the theory of non-classical orthogonal polynomials. In Section III 
we describe the general properties of symmetric a-Ensemble revealing the fact that a = 1 is a special point associated 
with peak formation in the density of states. This special case of border level confinement is investigated analytically 
in the Section IV in the framework of the theory of orthogonal polynomials, exploiting the properties of the symmetric 
Pollaczek polynomials. In Sections III and IV we find the density vi^j (x) of levels, two-point kernel Kn {x, y), locally 
rescaled two-level cluster function I2 (s, s') and smoothed connected correlations [i^n {x, ?/)]f,o„ between the density of 
eigenvalues. It will be shown that for strong and border level confinement the rescaled local two-level cluster functions 
as well as the smoothed connected correlations follow the universal forms which are typical of random-matrix ensembles 
with the Wigner-Dyson level statistics Section V contains conclusions. 



II. SYMMETRIC a-ENSEMBLE: CASE OF STRONG LEVEL CONFINEMENT (a > 1) 



The rigorous treatment of a-Ensemble with strong level confinement has been made possible by the recent devel- 
opment of the theory of non-classical orthogonal polynomials. The relevant polynomials P^"'' (x) are known as those 
orthogonal with respect to the Freud weights , . 

All the n-point correlation functions can be expressed through the two-point kernel 
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where fcjv is a leading coefficient of the orthonormal polynomial (x). The pointwise asymptotics of -P^"* (x) for 
large N have the form 
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Here 
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is the Nevai-Ullman density. 



D {N, a) = N^/"D (a) , D (a) = 
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and 



kN^i/kN = -D{N,a). 
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The straightforward calculations based on Eqs. (||) - (^) lead to the following expression: 
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(x)) cos (y)) - cos (x)) cos ($(,"1^ (y) 



Despite of this function is rather complicated, it can be rewritten in the universal form if one supposes that \x — y\ 



is much smaller than the scale <; of characteristic changes of the mean level density, <; ~ i^]^ 



Assuming also that both x and y stay away from the band edge D (iV, a) of the spectrum and making use of the 
asymptotic identity (x) — <i>j^^ (x) — arccos (x/D {N, a)), we obtain in the leading order in 1/iV 



where 
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sin [t^v^n - y)) 
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is the local density of levels. Note, that Eq. ( |l0|) is valid on the scale \x — y\ which is much larger than the mean 



level-spacing (^'^tv^) 



Correspondingly, in the large- A'^ limit the two-level cluster function 
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y=y(s') 



being rewritten in the terms of eigenvalues measured in the local level spacing s — xv^^^ and s' — yv'^^^ , locally follows 
the universal form 



Y2 (s, s') = 



SVC? [tt (s — s') 

k (.-.')]' 
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irrespective to the value a > 1. 

As is well-known |l], the level-spacing distribution function P (A) can be expressed through the eigenvalues {A (A)} 
of the Fredholm integral equation where ^/Yz — sin [tt (s — s')] /n {s — s') stands for the kernel 
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P(A) 
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(15) 



Since Y2 (s, s') for a-Ensemble coincides with that for Gaussian Ensemble, it inevitably leads to universal Wigner- 
Dyson statistics. 

The connected correlations between the density of eigenvalues for x ^ y are given by [p| 



= -K 
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and known to oscillate rapidly on the scale of the band width D {N, a) . The smoothed correlation function which is 
useful for the calculation of integral characteristic of spectrum can easily be determined. Bearing in mind Eq. (^ we 
obtain after averaging over rapid oscillations 
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Equation ( |17| ) proves that smoothed correlations in a-Ensemble with strong level confinement follow universal form 

0,0- 

Note that along with the proof of universality of eigenvalue correlations we have obtained an exact formula for the 
density of levels in a-Ensemble. Thus, Eqs. (0) and (pi]) yield following expression for a > 1: 
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that can be rewritten in two equivalent ways by means of Hypergeometric functions, 
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where z ~ x/ D {N, a). 

In the particular case a = 2 Eq. ( pO|) recovers the famous Wigner's semicircle, i/^-* (x) = ■n^^\/2N — x^ . 

It is worth pointing out that Eqs. (|19[) and (20) obtained within the framework of the theory of orthogonal polyno- 
mials exactly coincide with the density of states calculated for symmetric a-Ensemble in the mean-field approximation 
. This circumstance justifies the validity of the mean-field approach for a > 1. For pure linear confinement potential 
(a ~ 1) asymptotic formula Eq. (0) fails, signaling that point a = 1 is a special one. 



III. GENERAL PROPERTIES OF SYMMETRIC a-ENSEMBLE 



The special character of the point a = 1 may be understood appealin g to the recent mathematical literature on 



the theory of orthogonal polynomials with respect to the Freud weights |15|. Noting that the density of levels i^n 
for random-matrix ensemble with confinement potential V (x) is related to an inverse Christoffel function A^'^ (x) = 
J2i^=io^ Pf (^) foi' polynomials Pi {x) orthogonal with respect to the weight w {x) — exp {—2V {x)) as 

(x) = e-2^(-)A^i (x) , (21) 

we conclude that the crucial changes occur in the A^-dependence of the level density at the origin [p^: 



(iV,a)-< InA^, a = l . (22) 




Thus, at the point a = 1 the functional dependence of the density of states on the number of levels differs from that 
both for < a < 1 and a > 1. The ratio 

(23) 

with vijg being the background density of levels, demonstrates all the more dramatic behavior. The i>i,g can be 
estimated as N/2D {N, a), where D {N, a) is the band edge for symmetric a-Ensemble. The estimates of Christoffel 
functions |15| ] allow us to relate D {N, a) to the maximal zero xin of the corresponding A'^-th orthogonal polynomials, 
xiN ~ N^i^, so that viig ~ iV^~^/". Such a definition of v^g is relevant not only in the case of strong level confinement, 
a > 1, when density of levels in the spectrum bulk is smooth, but also in the case of weak level confinement although 
the bulk level density is no longer constant. In the latter case the density of levels in the spectrum bulk is very 
well approximated by l/|a;|^~" 0|. Since in the bulk of the spectrum — eD{N,a) with < e < 1 we obtain 
i^bg (N.a) e°'''^D{N,a)°'^^ ^ N'^-^/°'. This estimates is fully consistent with the definition i^bg = N/2D{N,a) 
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given above. Then we immediately obtain from Eq. (^2|) that ^ (a) = for < a < 1. The case a > 1 can be treated 
with the aid of Eq. ( po|) which yields the density of levels at the origin {x = 0) 

so that ^ (a) = ^ (1 — l/a) for a > 1. Finally, we obtain 

w ^ / 0, 0<a< 1 

The function f (a) is plotted in the Fig. 1. 

Equation d2q) implies that the transition point a = 1 corresponds to the formation of a sharp peak at the spectrum 
origin which also holds for < a < 1 (soft confinement), but is absent in the case of strong confinement potential, 
a > 1. 



IV. BORDER LEVEL CONFINEMENT (a = 1) 



A. Pollaczek random-matrix ensemble (PRME) 



The formalism developed in Section II cannot be directly applied to the case a — I, since the relevant pointwise 
asymptotic formula Eq. (^) fails. Nevertheless, this difficulty can be avoided by choosing V (x) in the form 

y(^) (x) = - - ln«;(^) (x) = - ^ In h + j + V^''^ (0) , (26) 

where w^^^ (x) is the weight function for Pollaczek polynomials. Their basic properties are collected in Appendix A. 

The behavior of confinement potential y^'*') (x)can easily be obtained from Eqs. ( p^ ) and (A2). In the vicinity of 
the origin x = 0, |x| <C A, expansion of the Eq. ( p6| ) yields the quadratic potential 

V'^^^ (x) « V^^^ (0) + ix^^^i) (A) , (27) 

where 4'^^) (A) = J^kLo ~^ '^)~^ ^ trigamma function. The long-range behavior of F'-'^-' (x) follows from the 
asymptotic representation of |r (A + ix)\ when \x\ — > oo (see Eq. (A2)), and turns out to be 

vWix)^^\x\-(^X-^yn\x\-Vi'^ (28) 

with V^^ = Aln2. 

Thus, the confinement potential in PRME exhibits such a long-range behavior which involves, in particular, the 
precisely linear growth at large if one puts A — 1/2. The calculations show (see Fig. 2) that only small discrepancy 
between linear confinement potential Vl (x) — tt |a;| /2 and V^^^/^) (x) takes place in the region < 1, which is 
negligible as compared with the band edge D {N, 1) — N. This circumstance allows us to correlate PRME with the 
transition point a = 1 of the a-Ensemble mentioned in the Introduction. The advantage of the proposed new unitary 
ensemble of random matrices is that it can be treated exactly in the terms of orthogonal Pollaczek polynomials. [For 
the sake of convenience, the derivation of asymptotic properties of these polynomials is entered in Appendix B] . 



B. Density of states 



The density of states can be calculated making use of the asyniptotics of Pollaczek polynomials found in Appendix 
B. Namely, the density of states for eigenvalues of PRME reads |l[ 

.^^^ (x) ^ e--<^'(^)-^ (Pj^l, (x) Aplf) (.) - P(^) (.) ^ . (29) 
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Bearing in mind Eqs. (A2), (A4), (A5) and different asymptotics for P^^^ (x) near the origin of the spectrum (Eq. 
(B14)) and in its bulk (Eq. (B9)), we easily obtain for finite A in the limit 1: 



(x) = - (In (2N) 

TT 



Re-^ (A + ix)) , \x\ < V2N, 



(30) 



Ztt 



, 1< Ixl < A^. 



(31) 



Here ^ (z) = (djdz) InF (z) is a digamma function 



Equations (|30|) and (|3l|) lead to the conclusion that the density of states for PRME does not depend on the 
parameter A in the bulk of spectrum, whereas at the origin this A-dependence holds. Such a behavior of the density 
of states is not a surprise and known for the generalized Gaussian and Laguerre ensembles Jl6| . 

ks N ^ oo the density of states at the origin tends asymptotically to the value v^^^ (0) = (In (2A^) — (A)V7r. We 
note that this result is in agreement with Vorig {N, 1), obtained for strictly linear potential Vl {x), see Ec^. (p2|). The 

non-asymptotic formula for (0) that is 



(0) 



3=0 



1 



' (0) 



can also be obtained. From recurrence Eq. (Al) it follows that 



P. 



(A) 



2,ii (0) = 0, P^'^ (0) = (-l)^ 
Therefore we arrive at the relationship 



r(j + A) 
r(A)r(j + i)- 



(32) 



(33) 



(0) 
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TT 



r(j + i)r(j + A) 
j-'^ r(j + i)r(j + A + i) 



(34) 



Here Eqs. (A2) and (A4) were used, and [m] stands for integer part of m. 

Note that the density of states given by Eq. (^l|), being extended onto the whole interval x € [— A^, A^] of the 
eigenvalues of PRME, has logarithmic singularity at the origin, but it still remains integrable and obeys normalization 
condition 



+ N 



N 



dxv^^^ {x) = N. 



(35) 



Moreover, Eq. (^Tj) exactly coincides with the density of states for confinement potential Vl {x) = ■k\x\/2 which 
can be found within the mean- field approximation j7|. This circumstance is a strong evidence that the mean- field 
approach, which was proved to be valid for V [x) ~ " with a > 1 (see Section H), is justified for weaker symmetric 
potentials up to linear- like (except for the region closed to the origin of the spectrum). 

The density of states is presented in the Fig. 3, where an excellent coincidence is observed between analytical 
asymptotic expressions Eqs. (^o|), ( ^l|) and the density of states, calculated from the Hypergeometric representation 
of PoUaczek polynomials (see Eqs. (A2) and (A6)): 



(X) 



\T(\ + ix)\ 



ra=0 



(2A), 



2Fi(-n,A-Hia;;2A;2) 



(36) 



We also computed the density of states for pure linear confinement potential Vl (a:) using the well-known in the 
theory of orthogonal polynomials matrix representation for Christoffcl function p7[ | , 



Vt^(x) = -e-2^^(") det 
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Ml 
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MAf-l 
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where 



Mfc 



/+CX5 
x''exp{-2VL ix))dx = 
-oo 



r(fc + i). 



(38) 



The resuhs plotted in Fig. 4 constitute additional justifications of the use of confinement potential l/^^/^) (x) instead 
of initial linear one, Vl {x). 



C. Two-point correlators 

The two-point kernel can be determined as 



(39) 



Here, again, different asymptotics must be used for the bulk of the spectrum and for its origin. 
In the bulk of the spectrum, 1 <C x < iV, and 1 y < N, Eq. (B9) gives 



sm (<(.)) cos (.-(.)) 
3in ($(^^) (y)) cos ($(^^) (x)) ^ X (40) 



1/4 



2 \ 



where 



{x)^^ + {N + A) arccos (-^) - nxiy^j^^ (x) . (41) 



It can be seen that despite of this function has a rather complicated form, the local properties of the two-point 
kernel remain to be universal. Really, for \x — y\ that is much smaller than the scale of characteristic changes of the 

mean level density, ~ ly^j^^ (^dv'"^'^ / dxj ^ N, and both x and y at finite distance from the edge of the spectrum, 
the third term can be neglected. Then, one obtains 



K 



(A) 
N 



i^, y) = ^^^^ [^N^ (^) - (y)) • (42) 



Making use of the different representation for <I>^'' , 



{x) = J (1 + 27V) + A arccos (^) - ^ v^^^ (z) dz, (43) 

with i^^'' defined by Eq. (|3l|), we easily obtain the universal form of the two-point kernel in the limit N —> oo: 

{x, y) = ^ sin [ttv'^^^ [x - y)) , (44) 
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where v 



.(A) 

N 



i^^-* {{x + y) /2) is a local mean level density. Equation proves the local universality of the two-point 



kernel and is valid on the scale |x — y| which is much larger than the mean level-spacing 



In the vicinity of the spectrum center, \x\ <ti v'2N, and \y\ ^ v2N, the asymptotic formula Eq. (B14) and Eq. 
M) yield 



■n{y-x) 



sm 



(4^' (2/)- (x)) 



(45) 



{x) =x\n (27V) - arg T (A + ix) = tt 



/A) 



{z) dz, 



(46) 



where v^"^ is determined by Eq. (|3^). Since near the origin the scale <; of characteristic changes of the mean level 

density is of the order A, Eq. ( p5| ) can also be rewritten in the universal form Eq. (^) but with density of states 
corresponding to the center of the spectrum. Here, again, universal form of the two-point kernel is valid on the scale 

\x — y\ which is much larger than the mean level-spacing (^^iv^) ^ 1/lnA'^. 
Thus, we arrive at the conclusion that two-level cluster function for PRME, 
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(A)2 
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{x,y) 



y=y(s') 



being rewritten in the terms of the eigenvalues measured in the local level spacing s 
takes the universal form Eq. (|l^) for any finite A. 

The connected correlations between the density of eigenvalues are determined as 



Xl' 



(A) 
N 



and s' 



(47) 
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= -K 



(A)2 
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(48) 



a X y and oscillate rapidly on the scale of the band width D {N, 1) = TV. The smoothed correlation function can 
easily be determined. Bearing in mind Eq. (^0|) we obtain in the bulk of the spectrum after averaging over rapid 
oscillations 



{x,y) 



1 



D{N,lY-xy 



27r2 {x — y) 



(49) 



The smoothed correlations near the spectrum origin can be calculated by means of Eqs. 
to be 



1 



con 2Tr^ {x — y)' 



r, X y^y. 



|) and (14^), and turn out 



(50) 



Equation (50) is the limited case of the Eq. ( |49| ) when both and \y\ are much smaller than the band edge D (iV, 1) 
Thus, Eqs. (M) and (pG) prove that smoothed correlations in PRME follow universal form. 



V. CONCLUSION 



We have considered the general properties of symmetric a-Ensemble with confinement potential V (x) ^ \x\" 
and established that the phenomenon of the peak formation in the density of states takes place. Namely, we have 
demonstrated that this sharp peak occurs for < a < 1 and is absent in the opposite case a > 1. In this sense the 
point a = 1 is a transition point associated with border level confinement. It has been shown that this transition 
point may be explored by means of slight changes in strictly linear confinement potential Vl (x) near the spectrum 
origin leading to the potential (x) = — ^Inw*^^/^^ (x) — ilncosh(7rx) which is a special case of more general 

confinement potential Eq. ( p6| ) connected to the symmetric Pollaczek polynomials. 

We have calculated the density of states, two-point kernel, two-level cluster function and smoothed correlations 
of the density of eigenvalues in the large- limit for a-Ensemble with strong and border level confinement. It has 
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been shown that both properly rcscalcd two- level cluster function and smoothed correlations of level density take the 
universal forms which are typical of random-matrix ensembles with the Wigner-Dyson level statistics. 

We have also demonstrated that the mean-field approximation is valid for calculation of the level density in symmet- 
ric a-Ensemble with strong level confinement, a > 1. In the case of border level confinement, a = 1, the mean- field 
approach is proved to be justified in the bulk of spectrum, failing near its origin. 

We would like to stress that treatment presented in this paper is rigorous and does not appeal to commonly used 
conjectures and approximate methods. 
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Appendix A. Definitions and basic properties of symmetric Pollaczek polynomials 

Symmetric Pollaczek polynomials Pn^^ (x) are determined by the recurrence equation [ p8[ , p9[ | 

nPi^Hx)-2xP.j^\{x) + in-2 + 2X)pj^%{x) = 0, n= 1,2,3,... (Al) 

with P'i'^ix) = 0,pI^^\x) ^ 1, and A > 0. These polynomials are orthogonal in the interval ~oo < a: < oo with the 
weight function 

to \ (k + X) J 



SO that 



dxPi^'' {x)Pi^'> {x)w^^'> (x) = Snmhi^^ , (A3) 



^ r(n + 2A) 
T{n+1)- 

From recurrence equation the leading coefficient of Pi^'' can be found: 



The following Hypergeometric representation holds for Pollaczek polynomials: 

pW (x) = -^^^ 2F1 (-n, A + ix; 2A; 2) exp (i7rn/2) . 



n! 

Lastly, we present the generating function that reads 



fe=0 



F{x,w)^y Pt^\x)w'' ^ (^^] |u,|<l. (A7) 



+ iwj (l + u,2) 



(A4) 



2"" 

kn = 7- (A5) 

r(n + l) ^ ^ 



(A6) 
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Appendix B. Asymptotic formulae for Pollaczek polynomials 



1. Formula of Plancherel-Rotach type 

We start with Eq. (A7) which being reversed yields 

2m J 



(Bl) 



Here the integration is extended over a contour 70 that encloses the origin w = Q and does not intersect the branch 
cuts [i, +ioo[ and ]— ioo, —i] associated with the singularities of generating function. Choosing x = [n + l) cosO with 
e < < TT — e, where e is a fixed positive number smaller than n/2, we rewrite Eq. (Bl) as 



/-yn 1 



dw 



70 {i + w-^y 



■ exp 



(n + 1) I i cos 6' In ( '■ — ) — Inw 

' ^ 1 + iw 



For the calculation of this contour integral in the limit n ^ 1 the method of steepest descent 
The saddle-point condition is 



i cos 9 In 



1 — iw 
1 + iw 



0, 



whence 



.,±161 



(B2) 

can be applied. 

(B3) 

(B4) 



and the original contour 70 of integration enclosing the origin must be deformed to pass through the points w 
along the directions of steepest descent and to avoid intersections with the branch cuts. 

The contributions of both saddle points to the leading term of the asymptotic expansion of the integral (B2) are of 
the same order. Therefore 



where 



S (wsp) = — \cos6\ =F i ( cos 6 In 



tan 

2 4 



(B5) 



(B6) 



and 



(l + u;^' (2cos0)^' 



Making use of Eqs. (B5) - (B8) we obtain the following asymptotic formula: 



(n+ l)^"^exp(7r|a;|/2) 



{2\x\f-iV^ l-{x/{n+l)f ^'^ 



(B7) 



(B8) 



TT { X 

sm I — + (n + A) arccos [ 



■In 



l-^l-{xl{n + l)y 



l + Jl-(x/{n + l)Y 



(B9) 
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2. Vicinity of the origin: \x\ <Si \i^2n 



At the origin a; ~ the saddle-point approximation used in the previous section to calculate asymptotic value of 
the integral (Bl) does not work, and the Plancherel-Rotach type formula fails. In this interval the Darboux method 
i pTt turns out to be fruitful. 

To obtain some asymptotic expansion for a polynomial P„ in accordance with the Darboux Theorem we have 
to expand the corresponding generating function F (x, w) in the vicinities of its singularities e*'^'-' on the unit circle 
Iwl = 1 into series of the form 



F {x, w)=Yl ^™ (1 - we-^^'-y'^"^'" . (BIO) 



J ak+mbk 

Cm' (i - we 

m=0 

Then the expression 



m— k 



furnishes an asymptotic expansion for Pn{x). 

For PoUaczek polynomials the singularities of the generating function Eq. (A7) occur at the points w = ±i, in 
whose vicinities generating function can be expanded as 

Fix,w) = 2-'^-f2(~^^/") (^^-yil±^wr-'^'\ (B12) 



m=0 

Then the expression 



leads to the asymptotic formula for Pollaczek polynomials. 

The leading term results from m = in Eq. (B13), so that for finite A and n 3> 1 



P^'^ (-) - (i) l^l^l 1- (2n) - argr (A + - Jn) 



(B14) 



The next term (to = 1) of asymptotic expansion Eq. (B13) is of the order O (a;^ + A^) + (A — 1)^^ /2fi^ • This 
circumstance imposes the restrictions |A| <C V^n and |x| <C V2n for Eq. (B14). 
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Figure captions 



Fig. 1. Function ^ (a) demonstrating the sharp peak formation in the density of states at the origin of the spectrum 

when a < 1 

Fig. 2. Confinement potential (x) +V,^^^^ for PoUaczek random- matrix ensemble with A = 1/2 (dashed line) 

and linear potential Vl (x) = tt /2 (solid line). The only small discrepancy takes place in the narrow region |a;| < 1 

Fig. 3. Density of levels for PRME. Solid line: asymptotic formulae Eqs. (30) and (31). Dotted line: Hypergeo- 
metric representation Eq. (36). Parameters: A = 1/2, = 20 

Fig. 4. Density of levels for confinement potential Vl (x) = n \x\ /2 calculated from Eqs. (37) and (38) (dotted 
line) and for Pollaczek random- matrix ensemble with A = 1/2 (Eqs. (30) and (31), solid line). iV = 20 
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"^On leave from: Institute for Low Temperature Physics and Engineering, Kharkov 310164, Ukraine 
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